Validity of Wiedemann Franz law in small molecular wires 
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We report our investigations on the finite size effects of Lorenz number in a molecular wire. 
Using Landauer-Biittiker formalism, we find that for sufficiently long wires there are two validity 
regimes of Wiedemann-Franz (WF) law, the cotunneling and the sequential tunneling regimes, 
while in small systems only the first regime survives. We compare our results with standard Kubo 
formalism and explain its failure to detect the WF law in small systems. Furthermore, our studies 
on exponentially localized disordered wires show that Lorenz number value Lo predicted by WF 
law is obtained only in the cotunneling regime. Also, Lorenz number L exhibits typical distribution 
at different temperatures corresponding to different tunneling process. In particular, first order 
tunneling results in a low value of L whereas second order tunneling recovers the universal value Lq. 



PACS numbers: 73.23.-b, 73.63.Nm 

I. INTRODUCTION 

Thermoelectricity is attracting huge scientific interest 
on account of its apphcations in future energy produc- 
tion and utilizationir— . Much of the work has focussed 
on developing efficient thermoelectric materials to con- 
vert waste heat energy to electric current and in reverse 
to perform refrigeration. Compared to bulk materials, 
low dimensional systems have the potential to achieve 
improved thermoelectric efficiency owing to their highly 
peaked density of states and to the high density of inter- 
faces that could be used to reduce parasitic heat flow^"^'^. 
In this regard, studies related to thermal and electrical 
transport in nano scaled molecular wires have gained con- 
siderable attention. 

Thermoelectric conversion efficiency is characterized 
by the figure of merit ZT — TS'^a/n, where S is the 
thermopower and a and k are the electrical and ther- 
mal conductivity at temperature T. For practical appli- 
cations, it is required that ZT 2> 1. However, this is 
practically hard to achieve as the thermal and electrical 
conductivity are related by the Wiedemann-Franz (WF) 
law which states that the Lorenz number L — n/aT is 
constan1>ii. The constant value Lq = (7r^/3)(/sB/e)^ for 
non-interacting systems, where fc^ is the Boltzmann con- 
stant and e the charge carried by each electron. WF law 
follows from the single particle Fermi liquid (FL) theory 
which assumes that both electric and thermal current are 
carried by the same Fermi liquid particles at sufficiently 
low temperatures, so that the Sommerfeld expansion can 
be appliedi^. In a non-interacting system where this 
FL theory holds WF law is valid even in the presence 
of arbitrary disorder^'^. Studies in interacting systems 
showed that the law is violated largely due to the non 
FL behaviorii-i^. 

However, studies so far reported are limited to either 



one or two quantum molecules or in the thermodynamic 
limit. In mesoscopic physics, the thermodynamic limit is 
meaningless and one is interested in the transport prop- 
erties of finite systems. Finite size effects are expected 
to influence the properties considerably in a one dimen- 
sional quantum system. For instance, the thermopower 
S does not follow the relation S ccT as expected in the 
thermodynamic limiti^. Also, Lorenz number is shown 
to diverge in a finite size closed system^S. Above men- 
tioned works follow from Kubo formula, without explic- 
itly considering the connections to the baths. Recently, 
the transport properties obtained using the Kubo for- 
mula were shown to differ significantly from the Redfield 
Quantum Master Equation approach (QME) that takes 
into account the effect of bath'^^. In practice, it is re- 
quired that both ends of the molecular wire are in con- 
tact with baths, which are able to exchange charges and 
energy with the wire. The validity regimes of WF law in 
these finite size molecular wires is still an open question. 

Following the seminal work of Anderson'^^, disorder 
has played an important role in understanding the trans- 
port properties of mesoscopic systems. Moreover, disor- 
der has interesting effects on the properties of a finite 
system as the transport properties are affected consid- 
erably, depending upon its size. For instance, a disor- 
dered system shows ballistic transport if the localization 
length is very large compared to the system length^^. 
Conductivities of strongly localized systems decrease ex- 
ponentially with the system size?'*. However, it is not 
yet clear how the Lorenz number varies with disorder in 
finite systems. While considering disordered systems, it 
is worthwhile to study the statistical distribution of ob- 
servables. Studies in theses directions showed that the 
conductances of strongly disordered system follow log- 
normal distribution'^^. Hence, it is interesting to explore 
whether Lorenz number still holds this log-normal distri- 
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bution. 

In this paper, we investigate the vahdity regimes of 
WF law in a finite non-interacting molecular wire at- 
tached to reservoirs using Landauer-Biittiker formalism 
for phase-coherent quantum transport from one reservoir 
to the other. With this approach we obtain conductances 
rather than conductivities and the Lorenz number is re- 
defined as i = S/GT, with G and S electric and thermal 
conductances, respectively. We compare our results with 
two other commonly used formalisms, namely standard 
Kubo formalism and QME, in exploring the transport 
properties of finite sized systems. In the latter part of the 
work, the studies are extended to disordered systems. In 
particular, our results show that the validity regimes of 
WF law depend upon the finite size of the system even 
for a non-interacting system. For long clean wires, there 
are saturation plateaus of Lorenz number at the univer- 
sal value of (7r^/3)(fcB/e)^ at two temperature regimes, 
corresponding to the cotunneling and the sequential tun- 
neling regime. As the wire length is reduced, one of the 
plateaus vanishes. We also find that the standard Kubo 
formalism and QME approach fail, for any system size, 
to obtain WF law in the cotunneling regime. With disor- 
dered localized wires, WF law is valid only in the cotun- 
neling regime. In addition, Lorenz number shows typical 
distribution at different temperature regimes. 

The outline of the paper is as follows. We introduce the 
model system and discuss the Landauer-Biittiker formal- 
ism to calculate conductances and the Lorenz number in 
Secini Numerical results are discussed for a clean model 
in Sec. IIIII and for a disordered one in Sec. IIVI Finally, 
we draw our conclusions in Sec. [V] 



II. MODEL AND FORMALISM 

We consider a molecular junction formed by connecting 
a one dimensional molecular wire between two electrodes. 
The Hamiltonian of our molecular system is 

H = Hw + He + Hy/E'i 

Hw = - ^ + h.c), 

1=1 

He = ^ ^Ekjdl^dkj, 

j=L,R k 

HwE = '^{tkLc\dkL + tkRcl^dkR + h.c). (1) 
k 

The first term Hw corresponds to a wire of N sites with 
nearest neighbor interactions, the second term He to 
the two electrodes left (L) and right (R) and the last 
term H\ye to the wire-electrode couphng. The operators 
cl{ci) and dl.j{dkj) are creation (annihilation) operators 
of electrons in wire and electrode j, respectively. E^j is 
the energy of the kth electron in the jth electrode and 
tkj is its tunneling amplitude. Here, the electrodes are 



reservoirs of non-interacting electrons in equilibrium at 
some temperature T and electrochemical potential fj,. 

Below we outline the calculations of thermoelectric 
properties in our model using Landauer-Biittiker formal- 
ism. The current through the wire is due to the electrons 
tunneling from one electrode to another. An electron at a 
given energy E scatters at the junction and can be trans- 
mitted through it or reflected back. The probability to 
tunnel across the junction is given by the transmission 
coefficient T{E). Hence, the electric (Je) and thermal 
(Jg) current (from left to right reservoir) in the molecu- 
lar wire are given hy^ 

Je = lj dET{E)[}L{E)- fR{E)], 

J.^IJ dE{E - fi)T{E)[fL{E) ~ fniE)]. (2) 

Here, e is the electronic charge, h is the Planck's con- 
stant, fhiE) and fniE) are the electron Fermi distri- 
butions in the left (L) and right (R) electrodes with 
temperature T^^n and electrochemical potential ^il,r- 

UlAE) = {eM{E ~ ^iLA/kBTLA + ^y^ where ks is 
the Boltzmann constant). In this work, we are interested 
in the linear response of the system and hence assume 
that the differences A/x = — i^R and AT = Tl — Tr 
are infinitesimally small. Hence, in Eq. ([2|), the electro- 
chemical potential fj, fi]^ Ki and the temperature 
T^Tl- Tr. 

Using the nonequilibrium Green's function technique, 
the transmission coefficient can be expressed as, 

T{E) = Tr{TL{E)G\{E)TR{E)G,{E)), (3) 

where Tr is the trace, Gs{,E) = {E - Hw -^l-^r)^^ is 
the retarded single particle Green's function operator and 
Tl^r{E) = i[EL,_R(£;)-S[_^(.B)] are the level broadening 
functions. Ei and E/j are the retarded self energies of 
the left and right electrodes respectively. 

We assume wide band limit of the electrodes. Hence 
the level widths are energy independent and are given 
by 7j = 27r^j. \tkj\'^S{E — Ekj). Furthermore, we take 
IL = 1R = 7- Thus, Fl = ^c\ci, Tr = ^c\cN and the 
Eq. ([3]) can be rewritten as 

T(£;)=7'|(l|G.|A^)p. (4) 

Note that 7 is the coupling strength which physically 
measures the rate at which the electrons tunnel across 
the junction. 

Using the Taylor expansion, 

fdE) « fR{E) -f ^^Aa^ + ^AT, (5) 
in Eq. ^ , the response of the system is given by 
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where the Onsager coefficients Ln, L12, L21, L22 are 
given by 



ill 
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(7) 



The (isothermal) conductance G, defined as the elec- 
tric current under the application of the voltage A/z/e 
with no temperature gradient, is 



AT=0 



in 



(8) 



The thermal conductance ^, the heat current per unit 
temperature gradient for zero electric current, is 



and the Lorenz number L is 



L11L22 — L21L1 



(9) 




FIG. 1. Dependence of Lorenz number on temperature T 
for a wire of length A'^ = 110 coupled to the electrodes with 
strength 7 = 10"''. Here Lo = 7r^/3. Note that WF law is 
satisfied in two regimes at low temperatures. 



GT' 



(10) 



For a smooth function T{E), the Sommerfeld expansion^^ 
of the integrals in ([7]) to lowest order in ksT/Ep, with 
Ep Fermi energy, leads to the Wiedemann Franz law 



i = in = ^ — 



(11) 



In the following sections we investigate in detail the 
dependence of the Lorenz number L on the temperature 
T and coupling strength 7 for a clean and disordered 
molecular wire. 



III. CLEAN WIRE 

The transport mechanisms in the molecular system can 
be understood from the transmission function T{E) of 
the molecular wire. For a wire of N sites, there are 
N quantum states with discrete energies. The density 
of states and T{E) in the limit 7 — >■ consist of se- 
ries of delta functions corresponding to these energies. 
While coupling to electrodes, electrons can enter or leave 
the wire and hence these delta peaks are broadened due 
to the finite life time of the electrons. If the coupling 
is very weak then the density of states remain as delta 
peaks broadened by a factor proportional to the coupling 
strength 7 whereas for a strong coupling all the peaks 
merge. Depending upon its energy relative to the energy 
spacing between the different levels in the wire {^E), an 
electron can tunnel the junction mainly in three different 



ways: 1) Sequential tunneling: At temperature T 3> Ai5, 
the energy of the electron is very high and hence can tun- 
nel across the junction in a sequential manner by spend- 
ing a finite life time within the wire. Hence, the current 
in this regime is proportional to the coupling 7. 2) Co- 
tunneling: This is a second order tunneling process oc- 
curring at temperatures T <^ AE. In this regime, the 
current through the wire varies quadratically with the 
coupling strength 7. 3) Resonant tunneling: This occurs 
when the energy of an electron matches exactly one of 
the discrete energy levels in the wire. Under these con- 
ditions, the electron is transmitted with unit probability 
and the current through the system increases sharply. 

Transmission coefficient T(E) is a smooth function in 
the sequential and cotunneling regimes. Hence, the Som- 
merfeld expansion can be applied and we expect to re- 
cover WF law in these regimes. To investigate the valid- 
ity of WF law, we have plotted the variation of Lorenz 
number with temperature T for a wire of N ~ 110 sites 
in Fig. [TJ Here we present the results for fixed cou- 
pling strength 7 = lO"'*. (Throughout this paper we 
take h ~ e ^ k-B ^ 1- Also, we set = 0.) As shown in 
Fig. [1] we found two plateaus where the law is satisfied 
exactly i.e., when the ratio L/Lq is one. Our numerical 
analysis shows that these two valid regimes correspond 
to sequential and cotunneling process. Between these 
two plateaus, there is a region bounded by the temper- 
atures Ti and T2 where, lowering T, the Lorenz number 
increases initially due to resonance tunneling and there- 
after decreases when the mixing of higher order tunneling 
occurs. 
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FIG. 2. Variation of characteristics temperatures Ti and T2 
with the length A'^ of a molecular wire attached to the elec- 
trodes with coupling strength 7 — 10"''. Here Ti is the high- 
est temperature in the cotunneling regime and T2 is the low- 
est temperature in the sequential tunneling regime where the 
Lorenz number L=Lq. Between Ti and T2, WF law is not 
valid even at low temperatures. Note that a variation of the 
temperatures oc 1/A'' is obtained. 



In order to clarify our statement that WF l&w is satis- 
fied only in sequential and cotunneling regime, we study 
the variation of Ti and T2 with the number of sites N . 
Note that these regimes correspond to cases where the 
energy level spacing IS.E ^ T and Ai? ^ T respectively. 
The mean energy level spacing ls.E in the molecular wire 
decreases with the number of sites TV as Thus de- 

pendences of temperatures Ti,T2 oc 1/N are expected. 
This is indeed what we obtain in Fig. [21 Here, the top 
panel represents the highest temperature Ti in the cotun- 
neling regime where the WF law is satisfied. The bottom 
panel represents the lowest temperature T2 in the sequen- 
tial tunneling regime. Decreasing the temperature below 
T2 results in resonant tunneling and deviation of L from 
its constant value of Lq. The temperatures Ti and T2 
are calculated such that the ratio L/ Lq is one up to third 
decimal point. 

From Fig. [21 it is clear that the temperature T2 is 
0.075 for a chain of 40 sites. The decrease T2 oc 
suggests that T2 for = 10 should be ~ 0.18. In Fig. 
[3l we have plotted the Lorenz number ratio as a function 
of temperature for A'^ = 10 sites. From the figure, it is 
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FIG. 3. Same as in Fig. [Tjbut for a wire of length iV = 10. 
Note that here only in cotunneling regime WF law is satisfied. 



clear that at T = 0.18 the ratio £/Lo is 0.953. Also, 
Lorenz number increases linearly with decreasing tem- 
perature T without any saturation in the sequential tun- 
neling regime. Indeed, numerical results showed that the 
plateau of constant Lorenz number seen in the sequential 
tunneling regime decreases with decreasing N and is al- 
most absent below N — 40. This is justified as it follows 
from the analytical derivation of WF law that the Lorenz 
number is obtained only at low temperatures T <^ Ep 
where Sommerfeld expansion holds. Note that the con- 
duction band width for our model is [—2,2]. Thus, the 
Fermi energy, the energy of the highest occupied energy 
level, Ep is of the order of 2 (see for instance ref.^^). 
Another important observation is that the magnitude of 
violation of the WF law is almost independent of length 
of the wire. 

So far we have discussed only the length dependence 
of Lorenz number ratio for weak couplings to the elec- 
trodes. However, the strength of the coupling also plays 
an important role in modifying the transmission spec- 
trum. Hence we have investigated the Lorenz number as 
a function of the coupling strength 7 in a wire of length 

= 100 in Fig. [4] As seen from the figure, the large 
variation of Lorenz number is smoothened for stronger 
coupling. This is due to the fact that a strong coupling 
broadens the transmission function such that different 
resonance peaks overlap. Hence, it is difficult to observe 
the sharp increase in current due to the delta peaked 
transmission function. 

Our numerical analysis shows that the temperature Tq 
at which L is maximum and below which mixing regime 
(mixing of first order and second order tunneling process) 
occurs linearly increases with 7. This is illustrated in 
Fig. [5] where Tq is plotted versus 7 for a wire of 10 sites. 
Thus, for a wire with 7 — > 0, the temperature Tq 0. 
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FIG. 4. Dependence of the Lorenz number on the electrode 
temperature T for a wire of length A'' = 100 with different 
coupling strengths 7. The straight line corresponds to the 
results obtained analytically for 7^-0. 



FIG. 6. Comparison of Landauer and Kubo formula for 
Lorenz number of a wire of length A'^ = 200. The wire is at- 
tached to the electrodes with a strength 7 = 10~*. Note that 
the results of Kubo formula coincide with that of Landauer 
formula only above a temperature Tk- 
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FIG. 5. Temperature Tq as a function of the coupling 
strength 7 for a wire of length A'^ = 10. Here Tq is the 
temperature at which the Lorenz number takes its maximum 
value. 



Our analytical calculations show that the Lorenz number 
diverges as for 7 — > (Details of the calculations are 
given in the Appendix). This perturbative (in 7) result 
is plotted as a straight line in Fig. S) Also, Tq decreases 
oc with increasing length of the wire. 

Another approach commonly used in investigating the 
transport properties of finite systems is the Kubo formal- 
ism. Here, details regarding the baths (in our case elec- 
trodes) and coupling to the system (wire) are neglected 



and only the steady state distribution 
used. These assumptions are justified 
the linear response of an infinite system, 
las for finite systems are derived by extr 
of iV -> 00. The electrical and thermal 
and K, can be calculated using Eqs. dH) 
a and k instead of the conductances G 
the Onsager coefficients given by 



of the system is 
for investigating 
However, formu- 
apolating results 
conductivities, tr 
and dH) but with 
and S, and with 



Lii = e^T[Dii5{uj) + crii(a;)], 
L12 = eT[Di25{uj) + 0-12(0;)], 
L22 = T[D225{lo) + 0-22(0;)]. 



(12) 



Here, 



ZN 



E 

i,k 



[i\ji\k){k\j.m\i), 



'^'™(^) = ^7Tt -. — 2^ 



-PEi 



ZN 



Ei^Ek 



x{i\jl\k){k\jr,,\i)S{iO-AE). (13) 

Ei and \i) are the ith Eigenenergy and Eigenstate of the 
system, AE = Ei — Ek, is the frequency, /3 = l/fc^T, 
Z is the partition function and ji and j2 are the charge 
and heat currents. Currents ji and j2 are calculated as 
ji = Je and j2 = Jq-fiJe where Je[q] = iY,fjl^[hi-i,di]. 
di is the number of electrons in the wire and ^2 = ^/-i is 
the local system Hamiltonian. Note that Hw — X]f ~^ 
where hi = cjci+i + h.c. 

Fig. [6] shows the ratio of Lorenz number calculated us- 
ing the above formula for our model with N = 200 sites. 
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FIG. 7. Variation of the temperature Tk with the length 
of the wire A'^. Here Tk is the temperature at which Kubo 
formula deviates from that of Landauer's. Note that Tk varies 
oc 1/N with the length of the wire A*'. 



For comparison, results obtained using the Landauer- 
Biittiker formalism are also plotted. The results were 
the same till temperature T = 0.03 after which the ratio 
decreases by using the Kubo formula. Let Tk be the tem- 
perature at which deviation occurs. We found that the 
temperature Tk decreases cc l/N with the increase in the 
length N of the wire. This can be understood as follows. 
Eqs. (|12p are exact only for infinite systems for which the 
partition function Z = J2i e~^^' 3> 1. For small molec- 
ular wires, Z can be large only for high temperature T. 
Fig. [7] shows the plot of the divergence temperature Tk 
with the variation of length of the wire N. It is clear 
from the Figs. [5] and [7] that the temperature T2 < Tk- 
This implies that Kubo formula only gives the results for 
sequential tunneling regime and hence can predict only 
one regime of WF law even in arbitrarily long wires. 

Finally, we have compared the results obtained by 
means of the Landauer-Biittiker approach with the Red- 
field QME2ii21i^. By construction QME is first-order 
perturbative in the coupling 7 and hence reproduces the 
results in the perturbative regime of the Landauer for- 
mula. This regime is bounded by the temperature Tq 
from below. Since Ti < Tq, QME cannot reproduce WF 
in the cotunneling regime. 



IV. DISORDERED WIRE 

In this section, we discuss the validity of WF law in a 
disordered wire. We model the disorder by introducing 
on-site energies with randomness. The Hamiltonian of 
such a wire is 



Hdw = - ^ + h.c.) ^ tic\ci, (14) 

i=l 1=1 

where are random numbers uniformly distributed in 
the interval [— M^, VF]. In one dimension, even for arbi- 
trary small disorder strength, the system becomes ex- 
ponentially localized and exhibits insulating behavior—. 
Conductances of disordered wires decrease exponentially 
with the length of the wire as G = Goe~^^^, S — 
Spe"^/'^, where f is the localization length. Also, dis- 
tributions of the conductances are log-normal parameter- 
ized solely by its mean value. However, these conclusions 
are true only at T = K. At any non-zero temperature, 
the exponential decrease is not apparent as the electrons 
can hop from one localized state to another—. 
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FIG. 8. Logarithmic average of Lorenz number L as a function 
of electrode temperature T for a wire of length = 100 with 
disorder strength W — 5 and coupling strength 7 — 10~*. 
The straight line corresponds to the case of clean wire. Here 
regime (1), (2), (3) and (4) corresponds to cotunneling, mix- 
ing, resonant tunneling and sequential tunneling respectively. 
Note that only in the cotunneling regime the results are the 
same for both wires. 

Localization length is related to disorder strength as 
^ ~ IhjW^^^ . We restrict our analysis to wires with 
^ ^ iV so that the system is insulating at all tempera- 
tures. For this we consider a wire of length — 100 with 
disorder strength W = h. Localization length of the wire 
is ^ ~ 1. As the conductivities exhibit giant fluctuations 
for different samples, we take logarithmic averages. De- 
pendence of logarithm of Lorenz number (InL) on tem- 
perature T for this model is depicted in Fig. [51 The 
values are obtained by taking average over 1500 disorder 
realizations. Logarithm of Lorenz number = 7r^/3 
is equal to 1.1908 (in units where e ~ kB = 1). From 
the figure, it is clear that the plateau of constant Lorenz 
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number Lq is recovered at low temperatures. This corre- 
sponds to cotunneling regime as indicated by (1) in the 
figure. As the temperature is increased, mixing regime 
[region (2) in Fig. |8] is reached where the Lorenz number 
decreases with increase in temperature. The decrease is 
apparent till the temperature where the resonant tunnel- 
ing occurs [regime (3)]. Further increasing temperature 
increases the Lorenz number until the sequential tunnel- 
ing regime. Finally, in the sequential tunneling regime 
as shown by (4) in the figure, Lorenz number decreases 
linearly with temperature. 

To better understand the differences in the variation 
of Lorenz number for the disordered and undisordered 
wire, we have plotted the variation of logarithm of Lorenz 
number L for clean wire as a straight line in Fig. |8l Our 
findings in this regard are summarized as follows: a) The 
Lorenz number for the disordered wire is always equal to 
or less than that of the clean wire at all temperatures, b) 
There is no saturation plateau of constant Lorenz number 
in the sequential tunneling regime. This is in contrast 
to the results in refi^ where it was pointed out that 
WF law is violated in the resonant tunneling regime and 
is valid far away from the resonance in the sequential 
tunneling regime. We note that due to the finite size 
of our wire sample to sample deviations from the WF 
law do not self average to zero in the higher temperature 
regime. Indeed, we numerically observe that the variance 
saturates to a constant value and is discussed in the latter 
part of this section, c) The temperature Ti at which 
the Lorenz number Lq is recovered is shifted to lower 
temperature. This follows from the argument that for 
^ ^ N, disorder reduces the width of the transmission 
peaks and the transmission function is smooth only at 
low temperatures compared to the clean case. Therefore, 
the Sommerfeld expansion substantiating WF law is valid 
only at lower temperatures than in the clean case. 

So far we have focused only on the mean value of the 
Lorenz number. However, in a disordered system with 
localization length much smaller than the system length 
i.e., ^ <^ N, fiuctuations can be as large as the average 
value. Under these conditions, only a statistical distribu- 
tion provides meaningful information about the system 
properties. It is well known that in these highly local- 
ized systems, conductance distribution is log-normal i.e., 
InG (or InS) follows normal distribution at T = K. In 
the latter part of this section, we analyze in detail the 
distribution followed by Lorenz number in the different 
tunneling regimes. 

1) Cotunneling regime: It is clear from Fig. [5] that 
in this regime we reproduce the WF law. From our nu- 
merical analysis we found that fluctuations of the con- 
ductances are almost twice as that of the average value. 
This has been advocated as the evidence for insulating 
behavior in one dimensional system'^^ and thus we ensure 
that disordered wire is indeed localized. We find that the 
fiuctuations of both conductances are almost perfectly 
correlated and hence the Lorenz number Lq with almost 
zero variance is obtained. This is clear from Fig. [3] (a) 
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FIG. 9. Distribution of Lorenz number L at different tem- 
perature regimes for a wire of lengtli A'' = 100 with disorder 
strength W = 5. Different panels correspond to temperatures 
a) T = 0.00001 in the cotunneling regime b) T = 0.001 in the 
mixing regime c) T = 0.01 in the resonant tunnehng regime 
d) r = 10 in the sequential tunneling regime. Note that in 
the panels (c) and (d), logarithm of Lorenz number is plotted. 



where the distribution at temperature T = 0.00001 is 
shown. A delta peak around the value 7r^/3 = 3.2898 is 
obtained at all temperatures in this regime. 

2) Mixing regime: In this regime we found that the 
logarithms of conductances still follow normal like dis- 
tribution with variance almost twice as that of average. 
However, the skewness is non-zero. Moreover, the fluctu- 
ations of both conductances are not perfectly correlated 
and therefore the constant value of Lorenz number Lq is 
not recovered. Interestingly, we obtain a bimodal distri- 
bution for the Lorenz number. This is shown in Fig. [9] 
(b) for temperature T — 0.001. One of the two peaks cor- 
responds to the value 7r^/3 while the other corresponds 
to a value near to zero. The peak around 7r^/3 is largely 
populated at low temperatures. By increasing the tem- 
perature, this peak is reduced while the peak with small 
Lorenz number is populated. This continues till the res- 
onant tunneling regime is reached where there is only 
one peak. Indeed, our results suggest that the two peaks 
correspond to first and second order tunneling process. 
For disordered wires, second order tunneling favors WF 
law with Lorenz number Lq whereas first order tunneling 
favors a small value of Lorenz number. 

3) Resonant tunneling regime: Here, similarly to the 
mixing regime case, the logarithm of conductances has a 
skew normal distribution (a normal distribution but with 
non zero skewness). In contrast to the mixing regime, the 
variances of the conductances decrease faster with tem- 
perature compared to their mean values and are almost 
of the order of the mean values. The distribution of the 
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FIG. 10. The top panel represents the distribution of the 
Lorenz number for a wire of length A'^ = 100 with disorder 
strengths = 5 in the sequential tunneling regime at tem- 
perature a) r = 20 and b) T = 1000 c). The bottom panel 
represents c) the variance (E^) and d) the skewness (ma/E'') 
of the Lorenz number distribution at different temperatures 
in the sequential tunneling regime. 



Lorenz number in this regime is always peaked near zero 
value with long tails towards the large value and can be 
approximated to a log-normal distribution. Hence, we 
plotted the distribution of the logarithm of L. A typical 
distribution in this regime at temperature T = 0.01 is 
shown in Fig. |9] (c) and is in accordance with our expec- 
tation. 

4) Sequential tunneling regime: Our results in this 
regime indicate that the logarithm of the conductances 
has a distribution similar to that of resonant tunneling 
regime. However, the variance initially decreases with in- 
creasing temperature and thereafter saturates. The dis- 
tribution of the Lorenz number is still peaked around very 
small value. Hence, similarly to the case of resonant tun- 
neling regime, we plotted the logarithmic distribution in 
Fig. El^d). Here, a skew normal distribution of increasing 
variance and decreasing skewness is obtained as the tem- 
perature is initially increased. However, after the initial 
change, both the variance and the skewness saturate to 
a constant value as the fluctuations of the conductances 
also saturate. Thus, a distribution invariant with tem- 
perature emerges and is demonstrated in the top panel of 
Fig. [TOl Here the panels (a) and (b) correspond to tem- 
perature T — 2Q and 1000 respectively. In order to bet- 
ter clarify the emergence of the invariant distribution, we 
have plotted the variance (E^) and the skewness (ma/E^, 
where is the third moment about the mean) of the 
distributions in the bottom panel. It is clear from the 
figure that the variance saturates around the value 15 
whereas the skewness fluctuates around —0.85 at high 
temperatures, indicating a temperature invariant distri- 
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V. CONCLUSIONS 

Using the Landauer-Biittiker formalism, we investi- 
gated the validity of WF law in finite sized molecular 
wires with and without disorder. For a clean system, we 
found that the validity regimes of WF law depend upon 
how an electron tunnels across the wire. In particular, 
Lorenz number Lq = (''''^/3)(fcs/e)^ is obtained in the co- 
tunneling and the sequential tunneling regime as long as 
the temperature in these regimes is much less than Fermi 
energy Ep. For wires with length TV < 40, the resonant 
tunneling occurs for temperatures larger than the Fermi 
energy Ep and hence WF law is valid only in cotunneling 
regime. We further compared our results with standard 
Kubo formula and with Redfield quantum master equa- 
tion and found that the two approaches diverge from the 
Landauer formula at particular temperatures Tk and Tq 
that decrease with the increase in the length of the wire 
N . The temperatures Tk and Tq are always larger than 
the temperature at which cotunneling occurs and hence 
it follows from our results that Kubo formula and QME 
will differ from the Landauer-Biittiker formalism even in 
the limit of infinite length of the wire. 

Furthermore, we explored an exponentially localized 
disordered wire using the Landauer-Biittiker formalism. 
Here, even for wires of length N = 100, WF law is valid 
only at very low temperatures corresponding to the co- 
tunneling regime. In contrast to an infinitely long wire, 
there is no self averaging of the sample to sample devia- 
tions from the WF law at higher temperatures in small 
wires. Moreover, Lorenz number shows typical distribu- 
tions at different temperatures corresponding to different 
tunneling processes. A delta distribution peaked around 
the value Lq is obtained in the cotunneling regime while 
a bimodal distribution is obtained in the regime where 
the mixing of first and second order tunneling processes 
occurs. Logarithm of Lorenz number shows skew nor- 
mal distribution in the resonant and sequential tunnel- 
ing regime. In particular, we find that a distribution with 
constant variance and skewness emerges in high temper- 
ature regime. We infer from our results that first order 
tunneling favors a small Lorenz number whereas second 
order tunneling favors the universal value Lq in a disor- 
dered wire. 

Finally, we point out that we have not addressed the 
effects of interaction between the electrons in our model. 
Commonly used approaches in investigating the trans- 
port properties of interacting finite sized systems are 
standard Kubo formalism and QME. We have discussed 
numerically in detail the failure of Kubo formalism in 
determining the WF law in finite sized systems. Further- 
more, it follows from our numerical analysis that stan- 
dard QME approach also fails. Standard QME is de- 
rived by taking second order perturbative expansion of 
tunneling amplitudes t^j and is only linear in coupling 
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strength 7. Hence, the currents calculated usmg this for- 
malism is always of first order in 7. However, in the 
cotunneling regime currents varies quadratically with 7. 
A fourth order perturbative expansion of tunneling am- 
plitudes in QME indeed explains this regime^. Thus, we 
note that to investigate the Lorenz number in finite sized 
strongly interacting systems the standard QME has to 
be extended to include terms of tfj . 
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Appendix A: ANALYTICAL DERIVATION 

In this section, we analytically derive the electric ( Je) 
and thermal (Jq) currents for a wire weakly coupled to 
the electrodes (7 —J' 0). 

From the definition of the single particle Green's func- 
tion operator Gg , it follows that Gs is the inverse of the 
matrix 



E-i^ -1 



-1 




E 



V 



E -1 
-1 E + i'-. 



J 



Thus, {l\Gs\N) is the (1, N) element of the inverse of the 
above matrix. 

Consider a wire with N=l coupled to both left and 
right reservoirs. Here, {l\Gs\N) = 1/{E — 17) and hence 



Je 



67 

2Trh 



dE 



1 



E^+r 



■ihiE) ~ fniE)] 



(A2) 



Apart from the constant factor, the first integral is of the 
form 



/ = 



/ dE^^fiE) ^ J dEg(E), 



(A3) 



which can be evaluated using residue theorem. The first 
part of the integrand has two poles namely E± — ±17 
whereas the second part has infinite number of poles 
namely Mastubara frequencies at En = l-i- + ^(2?^ + l)7r 
with n e Z and (3 — 1/T (we set /cs = 1). 

To evaluate the integral, we consider two distinct con- 
tours, one on the upper half plane (Ci) and the other 
on the lower half plane (C2). Each contour runs from 
— i? to i? on the real axis and then comes back following 
semicircle of radius R. Thus, the integration domain is 
the sum of the following two parts: 



R 

R JCr 



R 

-R JC' 



(A4) 



where Cji (C*^) denotes the upper (lower) semicircle. 
Since the integrand function satisfies the 'big circle 
lemma' (Jordan's lemma), the contribution of the inte- 
grals on the semicircles is zero when i? — !• cx) and thus 
only the contribution from the real axis integration sur- 
vives. Applying the residue theorem, we get 



/ = - 2ni[resg{E+) - Tesg{E^) 

00 —00 

+ ^Tesg{En) ' ^ Tesg{En)] 

n—0 n— — 1 



(A5) 



The factor 1/2 follows from the fact that integrating on 
two contours real axis integration is encountered twice. 
Now the first two terms in Eq. (lASP are 



res5(S+) = hm g{E){E - E+) 



E^E- 



= lim 



and 



f{E){E-E+) 
E^E+ {E-E+){E- E. 
f{E+) 

[E+~E.y 



resg(i?_)= liin g{E){E - E^) 



(A6) 



f{E-) 
{E+ - E- 



(Al) Thus, their differences is 



res g{E+) - res g{E_) = 



1 

iE+-E^) 

2Rc[/(i^+)] 

{E+-E-) 



-Re([/(*7)]. 

17 



(A7) 



[f{E+)-f{E^)] 



(A8) 



The residue corresponding to the nth Mastubara fre- 
quency En is 

1 1 

^^^9iEn} - {eMl3{E-^,}] + iy\E=E„ 

1 1 



E^ + -f^ 13 {ejip[P{E - + I) 

=-1 

1 1 



(3E^+j^' 
Here we point out that 

1 



(A9) 



E-n-i = M + ^(-2(n + 1) + l)iTr 

= fi+^{~2n~l)i7r^E*n. (AlO) 



Also, 



n=-l 



^ resg(En) ^ resg(i;_„,_i) {n -m - 1) 

m=0 

00 

= J2resg{E:). (All) 



n=0 
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Therefore, the difference between the residues is 



Tesg{En) - Tesg{E-n-i) = - 



/3 



1 



1 



^--2Inr 



1 



Substituting all these results in Eq. (|A2[) we get 



Je = g^|^Re[/L(z7)]-Re[Mz7)] 



1 



Im , — Fi- 



n=0 



1 



(A12) 



(A13) 



Following the same steps of = 1, we obtain 

Jp — zvri < = — 

2iTh 2 \ i-i{Ei -E2 + n){Ei - 

2Re[fL{E2 + i^)- fRiE2 + i^) 

ij{E2 - El + ij){E2 - El) 

In the limit of 7 — 0, the current Je is 

^ hiEi) - fRiEi) 



Je 



67 

T 

+ 



(El - E2Y 

fL{E2)-fR{E2) 



{E2 - EiY 
Similarly, for wire of length N we get 



N 



k=l 



3^k 



[Eu 



E,? 



E2 



(A17) 



(A18) 



(A19) 



where EnL,R. = IJ'L,r + I3]^\i'^n+ l)i7r with n S Z. In the 
limit of weak coupling i.e., 7 ^ T, the above equation 
reduces to 

Je ^ ^[fUE = 0) - ME = 0)]. (A14) 

Note that in the last formula E ^ should be substituted 
by the dot's energy E = eg in the case the single dot 
Hamitonian H\y = egc^c. 
For a wire with A = 2, 



(1|G.|A^) 



{E - Ei+i^){E - E2+q)' 



(A15) 



where Ej = 2 cos{[7r/(A^ + l)]j} (j = 1,2). Thus, the 
current Jg is given by 



67 

27r^ 



dE- 



fL{E)~fR{E) 



\{E-Ei+i^){E-E2+i^) 



(A16) 



for the electric current and 



J i^ {Eu-ti)[fL{Ek)~ fR{Ek)] 



(A20) 

for the thermal current, with Ej — 2cos{[7r/(A^ + l)]j} 
(j = l,...,Af). Using the Taylor expansion in Eq. ([5|), 
the difference of the Fermi functions reads. 



fL{Ek)-fR{Eu)^ 



1 



AT 



Am, 

T ^ 



4cosh2[MM 
AT 



A/i 



(A21) 



The above expressions are used in Eqs. ([5]), © and 
for calculating conductances G and 5 and the Lorenz 
number L in Fig. |4l It is clear from the above analytical 
calculations that conductances G oc 1/T and S oc l/T^ 
and thus the Lorentz number L — oc l/T^. 
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